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The statistics of a passive scalar randomly emitted from a point source is investigated 
analytically. Our attention has been focused on the two-point equal-time scalar correla- 
tion function. The latter is indeed easily related to the spectrum, a statistical indicator 
widely used both in experiments and in numerical simulations. The only source of inho- 
mogeneity/anisotropy is in the injection mechanism, the advecting velocity here being 
statistically homogeneous and isotropic. Our main results can be summarized as follows. 
1) For a very large velocity integral scale, a pure scaling behaviour in the distance be- 
tween the two points emerges only if their separation is much smaller than their distance 
from the point source. 2) The value we have found for the scaling exponent suggests the 
existence of a direct cascade, in spite of the fact that here the forcing integral scale is 
formally set to zero. 3) The combined effect of a finite inertial-range extension and of in- 
homogeneities causes the emergence of subleading anisotropic corrections to the leading 
isotropic term, that we have quantified and discussed. 



1. Introduction 

Turbulent fields whose statistical properties are invariant under translation and rota- 
tion in space are particularly interesting to theorists. They represent a formally simple 
setup allowing to neglect additional 'complications' introduced by the boundaries of the 
system or external driving mechanisms. However, in most situations of interest, such 
complications might play a crucial role in determining the statistical properties of the 
turbulent field. For the sake of example, this is the case in channel-flow turbulence, where 
statistical invariance under translation and rotation (i.e. homogeneity and isotropy) is 
only restricted to a small region arou nd the center of the channel, while it is totally lost 
close to the walls ( Toschi et al. 19991 ). 



In the last few years, important achievements have been obtained, both for the method of 
analysis and the level of comprehension, in relat ion to the effects on the sma ll-scale statis- 
tics of anisotropic large-scale contributions (see Biferale fc Procacciall2005 . for a review). 



The situation is much less clear for what concerns the role of inhomogcneities in determin- 
ing the statistics of small-scale turbulence. In particular, it is not clear if inhomogeneities, 
activated by boundaries and/or initial conditions, might disappear at small scales owing 
to cascade processes, which tend to eliminate a detailed memory of the large-scale dynam- 
ics. In terms of scaling invariance of small-scale statistics, it is unclear if the presence of 
inhomogeneities might lead to different scenarios with respect to the case where inhomo- 
geneities are absent. Understanding the above aspects has relevant consequences in appli- 
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cations related, for example, to small-scale subgrid parameterizations. Indeed, modern ap- 
proaches to closure problems commonly use scaling exponents as the basic ingredients to 
build subgrid scal e models. The best examples are the fractal dScotti fe Meneveau 1997 ) 
and multifractal ( Basu. Foufoula-Georgiou fe Porte- Age! 2004) interpolation schemes. 
How to define scaling exponents in the presence of anisotropies became clear only very 
recently (jBiferale fc Procaccial 120051 ). The situation is different in the presence of inho- 
mogeneities, where, up to now, also basic questions related to the existence of scaling 
behaviour, and thus of scaling exponents, still have to be seriously addressed. 
The main aim of our paper is to give quantitative answers to the above questions. To 
do that, we will focus our attention on passive scalar turbulence, where the scalar is 
randomly emitted from a point source. The injection mechanism is thus intimately in- 
homogeneous and the question on how this inhomogeneity eventually reflects on the 
small-scale scalar statistics can be addressed. In order to carry out the study in analyt- 
ical terms, we will assu me a white-in-time, ho mogeneous, random process to model the 
velocity-field statistics (jKraichnan I1968L Il994h . 

The paper is organized as follows. In §[2] we formulate the problem in the context of the 
Kraichnan advection model with inhomogeneous forcing and we adopt simple mathe- 
matical techniques to obtain the equation for the two-point equal-time scalar correlation 
function, whose general solution is provided in § [3] for a point-source emission. In § U 
we focus on the process of local cascade which represents an interesting example of 
persistence of inhomogeneity at small scales. Section [5] deals with the interplay between 
inhomogeneity, anisotropy and finite-size effects, which provides a correction to the domi- 
nant isotropic behaviour. Conclusions and possible future developments follow in § [6j The 
appendix is devoted to the reformulation of the problem in a finite box in the presence 
of periodic boundary conditions. 



2. Basic equations 

2.1. Kraichnan advection model 
Let us consider a passive scalar field 9(x,t) transported by a turbulent flow: 

d t + v ■ 89 = nd 2 9 + f , 



(2.1) 



where k is the molecular diffusivity. The incompressible velocity field v(x,t) is assumed 
statistically homogeneous and isotropic, whereas the source term f(x,t) is allowed not 
to be invariant under translations: a relevant example is provided by the emission of a 
tracer from a point source, located e.g. in the origin. 

Let us now specialize to the Kraichnan ensemble (jKraichnanl 1 196 1 Il994i where the 
velocity is a Gaussian, zero-average, white-in-time field with two-point correlation whose 
spatial behaviour is described by 



D tiu (xi,x 2 ) = £>o<W _ dvu>( x i ~ x z) ■ 
The second-order moment of the velocity increments is given by 

<W) = D ir e \(d + £ - 1)6^ - ^ 



(2.2) 



(2.3) 

for r = \r\ smaller than the integral scale of the velocity field (L v ), above which d^ y {r) 
saturates to an almost constant value whose order of magnitude is D\L^. Consequently, 
since the correlation D^ v (r) has to vanish for r — > oo, the relation Dq ~ D\L^ holds. 
Here, d is the space dimension (> 2) and £ is the scaling exponent, describing the degree 
of roughness present in the velocity field, lying in the interval (0, 2). 
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The two-point equal-time correlation function C(xi,x 2 ,t) = (9(xi,t)8(x2,t)} may be 
expressed as a function of the center of mass z = (x\ + X2)/2 and of the separation 
r = X\ — x 2 . In these coordinates the equation for the correl ation function C(r, z, t) 
follows from the ap plication of Gaussian integration by parts ()Furutsu 119631: Novikov 
19651 Donsker 1964 ). Its form reads 



d t C= [2K<J M „+d M „(r)] 



d 2 C 
dr^dry 



(D Q + 2 K )S^ + D^(r) d 2 C 
4 dz^dz u 



F 



(2.4) 



where F(r,z) represents the correlator (0(xi,t)f(x2,t) + 9(x 2 , t)f(x\, t)). 
Two cases appear to be quite relevant, also in connection with applications: the case of a 
constant emission from a point source and the one where the emission is random in time 
(but still punctual in space). The former case turns out to be quite cumbersome to be 
attacked by analytical methods and is still under investigation. Here, we shall focus on a 
Gaussian, zero-average, white-in-time forcing representing a random emission from the 
origin. Namely, f(x,t) = f (t)S(x) with (f(x 1 ,t 1 )f(x 2l t 2 )) = F Q S(x 1 )S(x 2 )S(t 1 -t 2 ), so 
that F(r,z) = Fo5(r)5(z). As we shall see in detail, this case is amenable to analytical 
treatment. 



2.2. Fourier transform and SO(d) decomposition 
Let us now come back to (|2.4I) . Fourier transforming it in z and defining 

C(r,q,t) =J d d ze- iqz C(r,z,t) , F(r,q) = J d d ze-' iqz F{r,z) , 

we obtain: 

d t C = [2k6^ v + d^ir)] ^ ^ q^C + F . (2.5) 

In the present case, the forcing-correlation transformed F = F 8(r) is independent of 
the wavenumber. 

Equation (|2.5p is differential only in r and is algebraic in the centre-of-mass wavenumber 
q. The second term on the right-hand side represents the inhomogeneous contribution 
and consistently vanishes for q — 0, which is equivalent to average all over the space. It 
is convenient to rewrite its r-dependent coefficient in the following way: 

_ (£>o + 2«)^ y + g^(r) 
4 

. (2.6) 

Substituting it back, it is clear that the last term in (|2.6p generates the only contribution 
in (|2.5p not invariant under rotations of r, because it gives rise to a scalar product 
between r and q that mixes different angular sectors. However, at separations r <C L v , 
a simplification is possible, since, in that case, the order of magnitude of d^(r) « 
is negligible with respect to Dq ~ D\L^. Therefore, D^ v (r) ~ DqS^ il , and the right-hand 
side of (|2.6p simplifies into —(-Do + k)(5 (LH/ /2. It is worth noticing that, when r is of the 
order of (or larger than) L v , a coupling between anisotropy and inhomogeneity takes 
place: we shall come back to this point in § [5j where the consequences of keeping the full 
form (|2.6|) into account will be discussed. Here, we concentrate on the sole case r <C L v 



Dq + k _ (d- + Q 
2 id 



5nu 
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in the stationary state with vanishing diffusivitjfj] and we can thus consider the simpler 
equation 



d MV {r) 



d 2 c 

dr^dry 



1 



-D q C + F = 



(2.7) 



A dimensional-analysis balance between the first and the second term in (|2.7p leads to 
the introduction of a new scale 

q 2 D l- 1/(2 -« 



_2(d- l)Di. 

which is associated to the strength of the scalar inhomogeneities and measures the sep- 
aration above which they b ecome relevant. A decomposition on the spherical harmonics 
(|Biferale fc Procaccia|[200i) . 

C(r, q) = VaJ2Cl,m(r, q)Yi, m ($) , F(r, q) = \/0^^ im (r, q)Y lim {$) , 



with $ denoting the solid angle associated with r and O its overall value, yields the 
following equation for C/ im (r, q) in each sector: 

r -( d -i )d y +£ -i drdi (d + Z-l)l(d-2 + l) r _ 2di (2 _ 6) £ + _ 

d — 1 H 
Note that, because of foliation on I and degeneration, we have dropped the dependence 
on the subscript m and we have introduced the rescaled forcing ipi (r) — Fi (r, q)/{d— l)Z?i 
(independent of q because of the punctual nature of the source) . 



3. General solution 

The general solution of (12.81) . as a function of r and £ q , reduces to the zero mode 
( Martins Afonso fc Sbragagli J2001I 



Ci(r; e q ) = w- Vo [AiK Ul (w) + B t I Vl (w)] 



(3.1) 



where w = 2(2 - O -1 ^/*,)* 2- ^ 3 and v x = [(d + t - 2) 2 + 4(d + i - l)l(d - 2 + l)/(d - 
l)]V*/(2-£). 

To determine the coefficients Ai and Bi , one can approximate the Dirac S by a Heavisidc 
0, exploiting the vanishing of ipi in all the anisotropic sectors I =/= 0: 



r -{d-l) r -(d-l) 

S(r) = S(r) = lim 6(L 



It must be stressed that, in this way, the point source is obtained as a limit of a forcing 



F or -( d -V 



■Q(L 



having positive Corrsin integral £ 


)n = [d d r F(r,q) = Fa. The investigation of the case 


Qn = (iFalkovich & Fouxon 12005 




Celani & Seminaral 20051 2006h will be left for future 



investigation. 

Studying the solution also for r < L, matching the solution C t and its first derivative in 
r = L , imposing regularity for small r and vanishing for large r (jMartins Afonso fc Sbragaglia 
2005), and eventually taking the limit L — ► 0, one finds Bi = VZ and 



A, = 



2-e 



Xd-S+2)/2 lim 



i / da,^(pK 0+1 ^,M = 5 Zi0 fc t -i^-d-S 



f Attention should, in principle, be paid to the limit of vanishing diffusive scale, but for the 
rough flows (£ 7^ 2) considered here no commutation problem arises with the limit of vanishing 
forcing correlation length L. 
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where u = w\ r=P = 2(2 - -1 OV4i) (2 ~ 0/2 . W = w lr=£ = 2 ( 2 - -1 (W) (2 ~ 0/2 and 
fcf = 2(2 - ^)-d/(2-€) /(d _ l)nr(^ + 1) (r(-) being Euler's function). 
In the pseudospectral space (r, q) the scalar-correlation transformed thus coincides with 
its isotropic projection and depends only on the moduli r and q (i.e. £ q ) as 

C(r;£ q ) = k^i-^-^r-^-^K^w) . 

Back to the physical space, the correlation is thus independent of the angle between r 
and z and is a function of r and z only: 



(J/ 2 



1 



n(2-Q 2 £>i 
AT 



47T 



z 2 r -(2-C) 



+ 1 



(3.2) 



hi = /fc t 2- 2 7T- 1 (2 - (4_e)/(2 -°r[2/(2 



£)]; for d = 3, k t = fc t 2- 1 7T- 3 / 2 (2 



For d = 2, /cj 

^(7-2?)/(2-?)r[3/2 + (i + 0/(2 - 0J- 

It is worth noticin g that the behaviour C ~ r -( d +t- 2 ) i w hich is typical of the homo- 
geneous situation l|Falkovich. Gawedzki fc Vergassolall200lf ) for r ^> L, is not observed 
in this case (even if L = 0), unless one integrates the correlation on the whole space, 
thus averaging out the inhomogeneity and defining the contribution in the homogeneous 
'sector'. In the pseudo-spectral space, this last operation is equivalent to consider q = 
(=> i q — > oo w = 0) and thus corresponds to keep into account only the leading term 
in the development of K Vo (w) for small arguments. 



4. Local cascade and small-scale persistence of inhomogeneity 

Recalling that the ratio D\ / Dq appearing in ()3.2() is of the order of , two opposite 
developments are meaningful, corresponding to small or large values of the quantity 

-vat ™ 

This defines a new characteristic length scale, 

C z = z 2/(a-€) L -f/(2-€) ; (4.2) 

dependent (monotonically) on z and whose meaning as the local integral scale of scalar 
turbulence will be clear shortly. 

For small s (i.e. r 3> C z ) the correlation is approximated by a power law in r, C ~ 
L v d ?/ 2 r -< i ( 4 -C)/ 2 + 2 -? , On the contrary, for large s (i.e. r <C C z ), a power law in z is 
found, and r appears only in subleading terms: 

C(r, z) ~ (9 2 (z)) - const, x e(z)r 2 ~ 6 , (4.3) 

Here, (9 2 (z)) = (Fo/Z3i)L« (d+5 " 2)/(2 ^ ) z 2 - d ( 4 -«)/( 2 -«) and e(z) = (9 2 (z)) / ' C 2 Z ^ repre- 
sent the scalar variance and the local dissipation, respectively. 

Expression (14. 3|) proves that, at scales smaller than C z (i.e. for r sufficiently smaller than 
z, according to the value of L v ), the increment of the two-point equal-time correlation 
function has the same behaviour shown by the homogeneous case (power law in r with 
exponent 2 — £: see e.g. Falkovich et al. 2001) at scales smaller than the forcing integral 
scale L. This result suggests that a cascade-like mechanism might thus be present also 
in this case, even if here r > L by construction. 

This possibility can actually be supported by the following considerations. The velocity 
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field sweeps the scalar, initially concentrated where it was released, and generates struc- 
tures which, Vxi, are correlated on the scale x\. This amounts to saying that correlations 
between each point X\ and the origin X2 = are created. In the centre-of-mass frame of 
reference, this means that in every point z (= Xi/2) a local cascade can then take place, 
starting from separations r sufficiently smaller than z. 

The quantity £ z thus plays the role of an effective local forcing correlation length: the 
word 'local' here refers to the fact that, while in the homogeneous situation the prefactors 
are expressed in terms of constant quantities, on the contrary in this case a dependence 
of (9 2 {z)) and e(z) on the point still persists. It is also worth noticing from (|4.2[) that, for 
very rough flows (£ — * 0), C z becomes proportional to z, consistently with the impossi- 
bility of a consistent definition of L v . On the contrary, for almost smooth flows (£ — ► 2), 
one finds C z ~ (z/L v ) 2 ^ 2 ~^ , critically dependent on whether the center-of-mass dis- 
tance from the source lies within the velocity correlation range. 

The aforementioned physical interpretation of local cascade can easily be supported 
mathematically by considering the physical-space counterpart of (|2.7p . i.e. the simpli- 
fied form of (|2.4p with the usual approximations r <C L v and n = = dt- 

d (r)^ + ^^ + F-0 (44) 

In the homogeneous case, the absence of any dependence on z gives rise to the convective- 
range balance 

d >- M £k - - F(r) ; (45) 



the derivative with respect to r of the left-hand side of (14. 5|) vanishes in the presence of 
a constant corresponding right-hand side, as is often the case. On the other hand, with 
point-source forcing, away from the origin, the balance (|4.4| has to be written as 

d (r)-^— - _£^_^£_ (46) 



but the vanishing of the derivative of the left-hand side of (|4.6| still takes place for r 
sufficiently smaller than z. This is shown in figure [TJ which moreover reflects how this 
interpretation has its validity limit affected by a change in the ratio r/L v appearing in 
the adimensional parameter s (|4.1[) . One should indeed remember that the three scales 
r, z and L v appear in a nontrivial way in s, whose magnitude is the key point for ap- 
proximation (14. 3|) and its consequences. 

The present case thus represents an interesting example of a situation in which inhomo- 
geneity persists at small scales, as the behaviour C ~ r 2- ^ observed at small r does not 
correspond to what was defined at the end of § |3] to be interpreted as the homogeneous 
counterpart (C ~ r -( d +£- 2 )y 



5. Finite-size effects and anisotropic contributions 

A comment is required about the relevance of the so-called finite-size effects. In other 
words, one would like to quantify the error deriving from the approximation r <C L v , 
which was used to simplify (|2.6p and thus to decouple inhomogeneity from anisotropy. 
This quantification becomes possible, if one proceeds in the following way. First of all, one 
should notice that, after the decomposition into spherical harmonics, no more foliation 
takes place. Namely, the equation for the isotropic sector is still a closed one (with the 
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appearance of a new term, according to the contribution within square braces in (|2.6p ). 



-Vd r (r d +^d r C )-l 



-(2-0 

q 



2d(d- l)D T 



and gives 



i + oi -r 



(5.1) 



but Cq now enters the equation for I — 2 as a forcing term (the I = 1 sector remains 
unforced because this procedure only couples even sectors, as can be deduced by de- 
composition (12. 6|) ). A simple power-counting operation is possible in Fourier space for 
r <C l q , where (specifying the order of the error in the development of K Vo (w) for small 
arguments) 

C oc i + (£) 2 ^ . (5.2) 

As a result, in this regime one easily obtains 



C 2 



r 



2-4 



C 



(5.3) 



Equation (|5.3[) shows that the first excited anisotropic sector carries a factor, with respect 
to the isotropic solution, given by the product between the corrections in (|5.1[) and in 
(|5.2p . Its interpretation is thus very simple and meaningful: at the lowest order, the most 
relevant anisotropic correction derives from the coupling of finite-size effects (0{r / L v )^) 
and of inhomogeneities (0(r/£ q ) 2 ~t). 

It can also be shown that, in the opposite situation (r ^> £ q ), Ci is still given by the 
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right-hand side of (|5.3j) but without the factor (r/£ q ) 2 ~^. This implies that, back to the 
physical space, the leading contribution is always the isotropic one, provided that r -C L v . 
The higher-/ anisotropic terms C4, Cq, etc. are indeed smaller and smaller, because they 
are forced by the (small) quantities C2, C4, etc. respectively. 

Such anisotropic corrections are expected to play a non-negligible role only when the 
scales r and z are comparable, but not when either is much greater than the other. An 
example of the former case is provided, for z — r/2, by the comparison between the 
situations z \\ r (where one of the two points in which the correlation is calculated lies 
on the source) and zlr (where both points are y2z away from the origin) : a difference 
must clearly exist, but cannot be caught by the isotropic function Co(r, z) and turns out 
to be subdominant. On the contrary, if r 3> z the two points are almost symmetric with 
respect to the origin, and if r -C z their relative separation is much smaller than their 
distance from the source: in both cases, a rotation of r with fixed z would change little. 
Of course, the problem is always invariant under rigid rotations of the whole space (and 
thus of both vectors r and z) around the origin. 



6. Conclusions and perspectives 

The dynamics of a passive scalar released from a point source has been investigated in 
this paper, as a prototype of inhomogeneity. Focusing on the Kraichnan advection model 
with Gaussian, white-in-time and zero-mean forcing, it has been possible to study analyt- 
ically the two-point equal-time scalar correlation function and to prove a persistence of 
inhomogeneity at small scales, in the spirit of the local cascade process described in § 0J 
Still to be understood in more detail is the interplay between inhomogeneity anisotropy 
and finite-size effects mentioned in § [5l 

An interesting open problem is the extension of our calculation to higher-order scalar 
correlation functions, with the aim of corroborating our results obtained for the two-point 
correlation function. 

Another natural extension of the present work would be represented by the study of a 
point source not satisfying the aforementioned hypotheses, e.g. a constant-in-time scalar 
emission. This situation is already under investigation, and a fully-analytical study should 
be completed by numerical results. A related issue is understanding how these results 
should change in the presence of smooth flows (£ = 2). 

Inhomogeneities were here limited to the forcing term, while the velocity field was as- 
sumed homogeneous and isotropic. It would be of interest to reformulate the problem 
with velocity ensembles not invariant under translation, so as to deal with a completely 
inhomogeneous situation. Moreover, whenever the Kraichnan model is used, the question 
of a possible extension to more realistic flows arises: numerical simulations would thus be 
required to test the present results with actual solutions of the Navier-Stokes equation. 
Lastly, a coarse-grained description might be applied to the present case, with the aim of 
finding large-scale closed equations and an explicit exact parameterization for the small 
scales, accounting for the effects played by inhomogeneities. The infrared limit of the 
theory for the passive sc alar has already been introduced with good results in the homo- 



geneous isotropic case ( Martins Afonso et all 120031: iMartins Afonso. Celani fc Mazzino 



l2004l ICelaniT Martins Afonso fc Mazzinoll2005l 120061 ). The extension to inhomogeneous 



cases is in order, in particular using the point-source problem as a paradigm. 
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Appendix. Periodicity and discrete spectrum 

It is interesting to reformulate the point-source problem in a finite d-dimensional box 
of side a with periodic boundary conditions, which is equivalent to consider an infinite 
d-dimensional grid of point sources with mesh size a. It can easily be shown that, upon 
Fourier transforming, the forcing-correlation spectrum is no longer continuum and flat, 
but is active only on a discrete set of wavenumbers (q^ = 2nk/a, Vfc £ *Z d ) with uniform 
intensities. In order to reconstruct the correlation function C(r,z) it is thus sufficient to 
analyse the discrete values C(r,qk) (we will only focus on the isotropic sector I = 0). 
For k = = qk, the pure homogeneous scaling behaviour C ~ r -( d +C~ 2 ) j s obviously 
found. For k ^ 0, such power law is replaced by Bessel functions, which can in turn be 
expanded in Taylor series, with a result analogous to (|5.2j) . That is, for each wavenumber 
the leading behaviour for small r is always given by the same homogeneous contribution, 
but for a fixed r the corrections to such term become more and more relevant with grow- 
ing modulus of gfc. When antitransforming, one would be tempted by extrapolating such 
leading behaviour from each mode and concluding that C ~ r -(<i+i-2) £ or sma vj r & \ so 
upon superposition. This is clearly not the the behaviour C ~ r 2 ^ is actually 

found for small r. In other words, it is not possible to exchange the Taylor and Fourier 
series, corresponding to the power expansion of the Bessel function and to the discrete 
antitransform respectively, because of the absence of uniform convergence. 
It is worth noticing that, in this consideration, we did not assume any power-law be- 
haviour for each mode (differently from what is usually done in each anisotropic sector 
upon SO(d) decomposition), coherently with the result that for each wavenumber we 
obtain a Bessel function, i.e. an infinite superposition of power laws. 
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